Electric Circuit Theory and the 
Operational Calculus 1 

By JOHN R. CARSON 
CHAPTER VI 

Propagation of Current and Voltage Along the 
Non-Inductive Cable 

THE principal practical applications of the operational calculus in 
electrotechnics are to the theory of the propagation of current 
and voltage along transmission systems. Of such transmission sys- 
tems the simplest is the non-inductive cable. The theory of the 
non-inductive cable is not only of great historic interest, relating as it 
does to Kelvin's early work on the possibility of transatlantic teleg- 
raphy, but is also of very considerable practical importance today, 
and serves as a basis for the theory of submarine telegraphy over long 
distances. We shall therefore consider the propagation phenomena 
in the non-inductive cable in some detail. 

The propagation phenomena in any type of transmission system 
are isolated and exhibited in the clearest possible manner when we 
confine attention to the infinitely long line, with voltage applied 
directly to the line terminals. Furthermore, as we shall see later, 
the solution for the infinitely long line is fundamental and can be 
extended to the more practical case of the finite line with terminal 
impedances. We therefore, in this chapter, shall confine our atten- 
tion to the case of the infinitely long cable with voltage applied directly 
to the cable terminals. 

Consider a cable of distributed resistance R and capacity C per 
unit length, extending from x = along the positive * axis. From a 
previous chapter (see equations (64) and (65) ), we are in possession 
of the operational equations of voltage and current; they are, for the 
infinitely long line, 

V= e-yl<rp y , (162) 

I= h^p e ~^ p Vo= ^ c i e ~^ p v °< (163) 

where a = x-RC, and V is the terminal cable voltage at x = 0. Let 
us now assume that the terminal voltage V is a "unit e.m.f."; then 

V=e-^ t (164) 

I = Rx^ / <*P e_V ^- ( 165 ) 

1 Continued from the October, 1925, issue. 
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The solution of (164) for V was considered in some detail in the 
preceding chapter; it is, by (129) 

V--J.I e -^=dr (166) 

where r = 4t/a = 4:t/x 2 RC. Series expansions of this solution were 
also given. Another equivalent form is, by (131) 






(167) 



This last form, recognizable also from inspection of the series expan- 
sion (132), is useful because the integral term is what is called the error 
function and has been completely computed and tabulated. 

Before discussing these formulas and the light they throw on propa- 
gation phenomena in the non-inductive cable, we shall derive the 
solution for the current. A very simple way of doing this is to make 
use of the differential equation (57) 

/— 4-5-7. 

Rdx 

Now from (166) and the relation 

9 _dr d 

dx dx dr 

8 v= 1 e~ l ' r d 4t 
dx -v/tt t \/t dx X °' RC 
2 P-i/r 



we get 



x\/t VV 
1= -4-7=-°-^.= J-4 e-"r. (168) 



whence 



It is worthwhile verifying the formula by direct solution from the 
operational equation (165). From formula (g) of the table of in- 
tegrals, we have 



fc = e -2Vx>VF^|! 
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Comparison with the operational equation shows that they are iden- 
tical, within a constant factor provided we put X = a/4. Conse- 
quently the solution of (165) is 



= ^-" /4 '=>li 



Rt e 



which agrees with (168). This, it may be remarked, is an excellent 
example of the utility of the table of integrals in solving operational 
equations. 

This formula is easily calculated for large values of / by expanding 
the exponential function; it is 

ifo v^L 1 \7) + 2!V7J -■••_]• 

The propagation phenomena of the non-inductive cable are there- 
fore determined by the pair of equations 

F 1 r^L *-l » f lHT e-"dr (109) 

\/ir >'o T\'T V 7T t/0 

and 

, 2 !l^lT rW (170) 

VttxR Vt \ ttRI 

4t 
where r = 4//a = -> D/ ^- 
X AC 

Now an important feature of these formulas is that the voltage 

at point re is a function only of -^, t; that is, of ±t divided by the total 

resistance and capacity of the cable from x = to x = x. The same 
statement holds for the form of the current wave: its magnitude, 
however, is inversely proportional to xR, or the total resistance of 
the cable up to point x. Consequently a single curve, with proper 
time scale serves to give the voltage wave at any point on the cable. 
Similarly a single curve, with proper time and amplitude scales, 
serves to depict the current wave at any distance from the cable 
terminals. These curves are given in Figs. 3 and 4. 

Referring to the curve depicting the current wave, we observe that 
it is finite for all values of *>0; consequently, in the ideal cable, the 
velocity of propagation is infinite. This is a consequence, of course, 
of the fact that the distributed inductance of the cable is neglected. 
Actually, of course, the velocity of propagation cannot exceed the 
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velocity of light. The error, however, in neglecting the inductance 
in the case of long cables is appreciable only near the head of the wave 
provided we confine attention to d.c. or low frequency voltages. This 
point will be discussed and explained more fully in connection with the 
transmission line. 

The current, while finite, is negligibly small until t reaches the 
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Fig. 3— Current in non-inductive cable (G=0) unit e.m.f. applied 
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Fig. 4 — -Voltage in non-inductive cable (£ = 0) unit e.m.f. applied 



value 0.2. In the neighborhood of this point it begins to build up 
rapidly; reaches at t = 2 its maximum value 



(0.429) 



yfixR V"2 VvxR 

and then begins to decrease, ultimately dying away in accordance 
with the formula 

2 1 f . 1 . 1/1\ 2 



"v/tt xR yV 



1-- 



1+hO 
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Its subsidence to its final zero value is very slow; for example, when 
t = 100 its value is still 

-2— (o.io). 

V ir xR 

Turning to the voltage curve, Fig. 4, we see that it is negligibly 
small until r reaches the value 0.25, at which point it begins to build 
up. Its maximum rate of building up occurs when r = 2/3, after 
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Fig. 5 — Power transmitted in non-inductive cable (G = 0) 

which it builds up more and more slowly. Its approach to its final 
steady value is in accordance with the formula 



V=l- 



V 7TT 



V 3x^2! 5r 2 7 



Even, therefore, when r is as great as 100, V differs sensibly from its 

ultimate value, unity, its value being 0.8876. 

x 2 RC 
Since the actual time is - r, it follows that the speed of building 

up is inversely proportional to the square of the length of the cable. 

The power curve VI is given in Fig. 5. V.I is the rate at which 
energy is being transmitted past the point x of the cable. 

The fact that the form of the current and voltage waves depends 
only on 4t/x 2 RC is at the basis of Kelvin's famous "KR" law, long 
applied to cable telegraphy and sometimes incorrectly applied to 
telephony. When the first transatlantic telegraph cable was under 
consideration, Kelvin attacked the problem of propagation along the 
non-inductive cable and arrived at formulas equivalent to (169) and 
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(170). From these formulas he announced the law that the "speed" 
of the cable, i.e., the number of signals transmissible per unit time, 
is inversely proportional to the product of the total capacity and total 
resistance of the cable (KR in the English notation). To see just 
what this means requires a little digression into the elementary 
theory of telegraph transmission. 

Telegraph signals are transmitted in code by means of "dots" 
and "dashes." The "dot" is the signal which results when a battery 
is impressed on the cable for a definite interval of time, after which 
the cable is short circuited. A "dash" is the same except that the 
time interval during which the battery is connected to the cable is 
increased. The "dots" and "dashes" are separated by intervals, 
called "spaces", during which the cable is short circuited. Now when 
the cable is short-circuited we may imagine a negative battery im- 
pressed on the cable in series with the original battery. Conse- 
quently the current in the cable, corresponding to a signal composed 
of a series of dots, dashes and spaces, will be represented by a series 
of the form 

I(t)-I{t-h)+I{t-h)-I(t-h)+I(t-U)- . . . (171) 

where, in the cable under consideration, /(/) is given by (168). h is 
the duration of the first impulse, t 2 — /i of the first space, h — h of the 
second impulse, etc. 
Now by (168) 

XR-V 7T V T XR-V 7T 

t is, of course, At/x-CR = ±t/KR (in the English notation). Now 
suppose that 

4/ 8 
T2 ~x 2 CR' 

Then the signal can be written as 



etc. 



-^^U(r)-<Kr-n) + <Kr-r 2 )- . . . } (172) 

Now if the relative time intervals n, t 2 . . . are kept constant (as the 
length of the cable is varied), the actual time intervals h, k . . . are 
proportional to x-CR or to KR, and the wave form of the total signal 
is independent of KR, when referred to the relative time scale t. 
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Hence, if T is the total time of the signal, T is proportional to x*CR 
(or to KR). That is to say, if the duration of the component dots, 
dashes, spaces of the signal are proportional to the "KR" of the 
cable, the wave form of the received signal, referred to the r time 
scale, is invariable, and the total time required to transmit the signal 
is proportional to the "KR" of the cable. Now the maximum theo- 
retical speed of transmission on the cable is limited by the require- 
ment that the received signal shall bear a recognizable likeness to 
the original system of dots and dashes: in other words there is a 
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Fig. 6 — Elementary telegraph signals in non-inductive cable 



maximum allowable departure in wave form between received and 
transmitted signals. If, therefore, the actual speeds of two cables 
are inversely proportional to their "KRs," the wave form will be the 
same. This establishes Kelvin's "KR" law. As a corollary, if the 
length of the cable is doubled the speed of signaling.is reduced to one- 
quarter, assuming the same definition of signals. 

The foregoing will be somewhat clearer, perhaps, if we refer to 
curves, 4, 5, 6, 7 of Fig. 6 which illustrate the distortion suffered by 
elementary dot signals in cable transmission. Curve 4 shows the 
dot signal produced by a unit battery applied to the cable terminals 
x 2 RC 



for a time interval t = 2 



'-, while curves 5, 6 and 7 are the cor- 



responding dot signals when the battery is applied for the time in- 

x 2 RC 1 x 2 RC 1 x 2 RC 

tervals : — ,-- . ~ * — : ; — and T '- — -. — . Any further decrease in the dura- 
4 2 4 4 4 

tion of the impressed dot, beyond that shown in curve 7, does not 
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affect the shape of the transmitted dot, which means that the cable 
speed has reached its theoretical maximum. These curves, it should 
be observed, can be interpreted in two ways. First, we can regard the 
length x of the cable as fixed and the duration of the impressed dot 
as varied. On the other hand, we can regard the actual duration of 
the impressed dot as constant and the length of the cable as varied. 
From the latter standpoint the curves illustrate the progressive dis- 
tortion of the signal as it is transmitted along the cable. 
The dot signal of relative duration T can be written as 

D = I(t), t<T 

= I(t)-I(t-T), t<T 

and the second expression can be expanded in a Taylor's series, giving 

If T is sufficiently short this becomes 

D = T I\t). (173) 

Hence when the dot signal is of sufficiently short relative duration 
T, the wave shape of the received signal is constant, I'{t), and its 
amplitude is proportional to the relative duration of the dot. 

This can be generalized for any type of transmission system : 
Let the dot signal be produced by an e.m.f. /(/) of actual duration T. 
Then the received dot signal, by formula (31), is 



D = j t £f(r)I(t-T)dr, t<T 

=W iT)I{t - T)dT > t>T - 

For /> T this becomes 

D = l\t) f r f(r)dr-I"(t) f T Tf(T)dT+ . . . 

Jo Jo 

and for sufficiently short duration T, we have approximately, 

D=I'{t) f T f(r)dr. (174) 

Jo 

Hence for a sufficiently short duration of the impressed e.m.f. the 
received dot signal is of constant wave form, independent of the shape 
of the impressed e.m.f., and its amplitude is proportional to the time 
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integral of the impressed e.m.f. These principles are of considerable 
practical importance in telegraphy. 

The leaky cable, that is, a cable with distributed leakage conductance 
G in addition to resistance R and capacity C, is of some interest. The 
differential equations of the problem are given in equations (70) ; the 
operational formulas for the case of voltage directly impressed on the 
terminals of the infinitely long line are 

y—g-xilcRp+RG y ot 



->lf 



f 

}f_ e -x\'cRp+Rc 7 o- 
R 



Writing CRx 2 = a and RGx 2 = (3, G/C = \, and assuming a "unit e.m.f." 
impressed on the cable, this becomes 

V=e~^«P+0, (175) 

*- \5v?+X e~V«Tw. (176) 

These equations are readily solved by means of the table and formulas 
given in a preceding chapter. 

But first let us attempt to solve the operational equation (175) for 

the voltage by Heaviside methods, guided by the solution of the 
operational equation 

V=e~^«p (124) 

ofthe preceding chapter. Expand the exponential function in (175) 
in the usual power series; it is 



/: 



(ap+B) ( a p+R)V<xp+B 



7=1- v«*+fl+£3^-'- w ' Yi + • ' < 177 ) 



Now discard the integral terms and write 

ap+& . (ap + PY . 



V=l 



1+ 



3! ' 5! 



Vap+P- (178) 



We have now to interpret the expression y/ap+fi. We have by 
ordinary algebra 



\ <*pj \ PJ (179) 
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Now identify y/p with 1/V \t in accordance with the Heaviside rule, 
and \/p with \dt. We get 



7=1- 



Now in the terms of the expansion (178) identify p" with d"/dt" 
and substitute (180); we get 

. 4(4+^(«^> 2 )+--. 

|— f x* (xo 2 (x/) 3 I (181) 

This series is hopelessly complicated to either interpret or compute. 
It is, in fact, an excellent illustration of the grave disadvantages under 
which many of Heaviside's series solutions labor. We shall there- 
fore attack the solution by aid of the theorems and formulas of a 
preceding section. The simplicity of the solutions which result is 
remarkable. 

The operational formula for the voltage is 

y =e -y/ip+0 m (175) 

Now the operational formula for the voltage in the non-leaky cable 
is (see equation (104) ) 

In order to distinguish between the two cases, let us denote the voltage 
in the latter case by 7°; thus 

p, e -VS#. (182) 

Now by theorem (VII) and equation (182) we have 



p+\ 

/>+X 
Now write (175) as 

V P p+X 6 

\ T p p+\ 



(183) 



(184) 
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It follows at once by comparison with (183) and the rule that \/p is 
to be replaced by ) dt, that 

V=(l+\rdt)v°e-\ (185) 

By a precisely similar procedure with the operational formula (176) 
for the current, we get 

I=(l+\ f l dt\l°e-*' (186) 

where 1° is the current in the non-leaky cable. Now by formulas 
(169) and (170) 



I^JJLe-*'*', (170) 

\ irRt 

which completes the formal solution of the problem. 

Formulas (185) and (186) are extremely interesting, first as showing 
the superiority of the definite integral to the series expansion — compare 
(185) with the series expansions (181) — and secondly as exhibiting 
clearly the effect of leakage on the propagated waves of current and 
voltage. We see that in both the current and voltage the effect of 
leakage is two-fold: first it attenuates the wave by the factor e~ ', 
(\ = G/C), and secondly it adds a component consisting of the pro- 
gressive integral of the attenuated wave. This, it may be remarked, 
is the general effect of leakage in all types of transmission systems. 
Its effect is, therefore, easily computed and interpreted. 

Formulas (185) and (186) are very easy to compute with the aid 
of a planimeter or integraph; or, failing these devices, by numerical 
integration. However, for large values of /, the character of the waves 
is more clearly exhibited if we make use of the identity 



dt 
Jo Jo Jl 

whence 



and 



f'dt = I dt- I 
Jo Jo Jl 

F=(l+X f°dt\ V°e- x <-\ I °° V°e~ x 'dt (187) 

1= (l + X i' X dt\l n e-*<-\ f *°I°e- h 'dt. (188) 
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The first two terms of these formulas are clearly the ultimate steady 
state values of the voltage and current waves, and can be deter- 
mined by evaluating the infinite integrals. A far simpler and more 
direct way, however, is to make use of the fact that the ultimate 
steady values of V and I are gotten from the operational formulas 
by setting p = 0. That this statement is true is easily seen if we 
reflect that the steady d.c. voltage and current are gotten from the 
original differential equations of the problem by assuming a steady 
state and setting d/dt = 0. 

From the operational formulas we get, therefore, 

(l+X /""<&) Ver*-e-1t-e-1™, (189) 

(l+\f"dt)l°e-*=^e~^ = ^e-^™. (190) 

Introducing these expressions into (187) and (188) respectively, we get 
F=e -W*c_x f°° V°e- Xl dt, (191) 

1= ^ e -WKi_x f~Pe-*dt. (192) 

The definite integrals can be expanded by partial integration; thus 
-X f a °V°e-*<dt= f Vde-*' 

J x d 

e~* 1 j Vdt. 

Continuing this process we get 

K=e -^_^( 1+ 4 + ^ + ..)^, (.93) 

Using the values of V° and 1°, as given by (169) and (170), it is ex- 
tremely easy to compute V and I, for large values of /, from (193) 
and (194). 

So far we have considered the current and voltage waves in re- 
sponse to a "unit e.m.f.," impressed on the cable at x = 0. It is of 
interest and importance to examine the waves due to sinusoidal 
e.mis., suddenly impressed on the cable, particularly in view of 
proposals to employ alternating currents in cable telegraphy. 
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We start with the fundamental formula 

x(t) = ^£f(t-r)h(r)dT 

= £f(t-r)h r (r)dr 

provided h{o) = 0, which is the case in the cable. 
If /(/) =sin ut, we write 

x 5 (t) =sin cat I cos wt.h'{t)dt 
Jo 

— cos o)t I sin ut.h'(t)dt. 

Similarly, if the impressed e.m.f. is cos ut, 

x c (t) =cos wt I cos u)t.h'(t)dt 

+sin ut J sin o)t.h'(t)dt. 



(194-a) 



(194-b) 



The investigation of the building-up of alternating currents and 
voltages, therefore, depends on the progressive integrals 

C = J o cos ut.h'(t)dt, 

„ (194-c) 

S= J sin ut.h'{t)dt. 

For the case of the voltage waves on the non-inductive, non-leaky 
cable these integrals, by aid of equations (169), become, if we write 
a/ = aco/4, 

r _ 1 C T e- x ' T cos o/r 

^ — ~~7 = I 7 = dr, 

(194-d) 



o 1 p g-'^sina/r 

\/x «^o t\/t~ ' ' ' 



where, as before, r=4t/a. 

For the current wave we have, by (170), 

2 fYi 1 \ e~ l,T cos gj?" 

vV^i U 2) V7 dr ' 



2 r r /l 1 \ g- 1/T sincoV 

S = y f /7xRJo \r 2) ryV dT - 



(194-e) 
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For small values of t and u these integrals can be numerically 
evaluated without great labor. Mechanical devices, such as the 
Coradi Harmonic Analyzer, are here of great assistance. In fact 
the Coradi Analyzer gives these progressive integrals automatically. 
It may be said, therefore, that a complete mathematical investiga- 
tion of the building-up of alternating current and voltage waves 
on the non-inductive cable presents no serious difficulties, although 
the labor of computation is necessarily considerable. One fact makes 
the complete investigations much less laborious than might be sup- 




30 x-4t/x*RC 



Fig. 7 — Non-inductive cable (G=0), building-up of alternating current. 

4 
Applied e.m.f. cos wt; w = 2ir ~Y57< 



posed. This is, if the foregoing integrals are calculated for a given 
value of co', the results apply to all lengths of cable and all actual 
frequencies u}/2ir, such that aw is a constant. Then if we double the 
length of the cable and quarter the frequency, the integrals are un- 
affected. 

The solid curve of Fig. 7 shows the building-up of the cable voltage 
in response to an e.m.f. cos «/, impressed at time / = (). The fre- 
quency u/2w is so chosen that w' = aw/4 = 2ir, and the curve is cal- 
culated from equations (194-b) and (194-e). The dotted curve shows 
the corresponding steady-state voltage on the cable; that is, the voltage 
which would exist if the e.m.f. cos w/ had been applied at a long time 
preceding / = 0. We observe that, for this frequency, the building-up 
is effectually accomplished in about one cycle, and that the transient 
distortion is only appreciable during the first half-cycle. 
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The case is very much different when a higher frequency is ap- 
plied. Fig. 8 shows the building-up of the alternating current in the 
cable when an e.m.f. sin cot is applied at time / = 0. The frequency 
is so chosen that o/ = a:co/4 = 107r. The outstanding features of this 
curve are that the initial current surge is very large compared with 
the final steady-state, and that the transient distortion is relatively 
very large. It is evident that the frequency here shown could not be 
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Fig. 8 — Non-inductive cable (G=0). Building-up of alternating current. 
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Applied e.m.f. sin ut; oj = 10tt 2Rr 



employed for signaling purposes. This curve has been computed 
from the steady-state formulas, and equations (160) and (161) for 
the transient distortion. 

If the applied frequency u/2ir is very high, the steady-state becomes 
negligibly small, and the complete current is obtained to a good 
approximation by taking the leading terms of (160) and (161). Thus 
if the applied e.m.f. is sin ut, and w is sufficiently large, the cable 
current is 

2 Id e~ l ' T 

by (1G0) and (170) while,, if the impressed e.m.f. is cos cot, it is 



2 /1\ 2 d 2 e-V* 



k wV 



by (161) and (170). Here w' = aco/4 and t=4//«. 
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CHAPTER VII 

The Propagation of Current and Voltage Along the 
Transmission Line 

We now take up the more important and difficult problem of 
investigating the propagation phenomena in the transmission line. 
The transmission line has distributed series resistance R and in- 
ductance L, and distributed shunt capacity C and leakage conductance 
G. It is the addition of the series inductance L which makes our 
problem more difficult and at the same time introduces the phenomena 
of true propagation with finite velocity, as distinguished from the 
diffusion phenomena of the cable problem. The cable theory serves 
very well for the problems of trans-oceanic telegraphy 8 but is quite 
inadequate in the problems of telephonic transmission. 

If J denotes the current and V the voltage at point x on the line, the 
well known differential equations of the problem are: — 

(4«)'"5* 



dx 
Replacing d/dt by p, these become 



(195) 



(Lp+R)I = -J-V, 



(Cp+G)V= -£-I. 



dx 

9 
d~x 



(196) 



From the second of these equations 

dv = i a 2 

dx " Cp+Gdx 2 
and substitution in the first gives 



/ 



{Lp+R){Cp+G)I = ^L (197) 



Similarly if we eliminate /, we get 

A-G\V= 

dx 



(Lp+R)(Cp+G) V=^V. (198) 



8 With the installation of the new submarine cable, continuously loaded with per- 
malloy, this statement must be modified. In this cable, the inductance plays a very 
important part, and is responsible for the greatly increased speed of signaling ob- 
tainable. 
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If we assume a solution of the form 

V=Ae-v x +Be* x 

where A and B are arbitrary constants, substitution shows that the 
solution satisfies the differential equation for V provided 

y°-=(Lp+R)(Cp+G). (199) 

From equation (19(3) it then follows that 

I= Lp\r Ue-y*-Bey*) 

(200) 

=$£±<i {Ae~y*-Bey*). 

7 

Now restricting attention to the infinitely long line extending along 
the positive x axis, with voltage V impressed directly on the line at 
x = 0, the reflected wave vanishes and we get 

v=v e-y\ 



i= c ±+G Vo e-y>, 

7 

Y- = (Lp+R){Cp+G). 



(201) 



Now let us write 



where 



7 2 = - 2 [(/>+p) 2 -* 2 ] (202) 



v = l/VLC, 
P ~2L + 2C 

__4__£ 

a 2L 2C 
Then setting V„ = l, the operational equations of the problem become 

7=J-7V<*+p)«-»», (203) 



v />y -v/^+p) 2 -^ 

Now consider the operational equation, defining a new variable F: 

F=P / =• (205) 

VO+p) 2 -* 2 
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It follows at once from our operational rules, and (203) and (204), that 

I = v(c+G f'dt^F, (206) 

V=s - V f'd£ d t. (207) 

Jo 8* 

Our problem is thus reduced to evaluating the function F, from the 
operational equation (205). This equation can be solved by aid of 
the operational rules and formulas already given. The process is 
rather complicated, and there is less chance of error if we deal instead 
with the integral equation of the problem 

-_= f °° F{t)e-»'dt. (208) 

-2 Jo 



e v 



ViP+py-o 2 



Now let us search through our table of definite integrals. We do 
not find this integral equation as it stands, but we do observe that 
formula (m) resembles it, and this resemblance suggests that formula 
(m) can be suitably transformed to give the solution of (208). We 
therefore start with the formula 






i- C™ e - p, JoWP-V)dL (m) 

1 Jx 



This, regarded as an integral equation, defines a function which is 
zero for t <X and has the value J (\/' 2 -X 2 ) for t^\, J being the Bessel 
function of order zero. We now transform (m) as follows : 

(1) Let \p = q and */X=/i. Substituting in (m) we get 

Vg 2 +X 2 J 1 

Now, in order to keep our original notation in p and /, replace q by p 
and /i by /; we get 

i_ = / e~i"J {Wt 2 -l)dt. (m.l) 

Vp 2 +x 2 J' 

(2) In (m.l) make the substitution p = q+n and then in the final 
expression replace q by p\ we get 
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x x 

(3) In (m.2) make the substitution p = — q and h= — t, and ulti- 

v v 

mately replace q by p and / 2 by / ; we get 

g-fV(P+/l 1 ) I +X 1 2 



3) 



where Xi= — X and ni = — n. (They are, of course, as yet, arbitrary 
x x 

parameters, except that they are restricted to positive values). 

(4) Now if we compare (m.3) with the integral equation (208) 

for F, we see that they are identical provided we get 

which is possible, since p>a. 

Introducing these relations, we have 

/ e-P'.e-P' I {Wt 2 -x 2 /v 2 )dt = . =. (m.4) 

Here I denotes the Bessel function of imaginary argument; thus 
Jo(iz)=I (z). 

It follows from (m.4) and the integral equation (208) that 

F(t) =0 for t<x/v, (209) 

= e ~i"Io(<rVP-*?/tP) for t = x/v. 

Having now solved for F=F(t), the current and voltage are gotten 
from equations (200) and (207). Thus 

1=0 for t<x/v, 

l~r r , (210) 

= J- F{t)+vG I F{t)dt for t = x/v. 
\ L Jx/v 

The corresponding voltage formula is 

7=0 for t<x/v\ 

n , , ax r i e-P T I 1 (<rVi 2 -x 2 /v 2 ) J . ^ . (211) 

- e ~px/v ^ / u . v '— '- dr for t = x v. 

V J x / V \/t 2 -X 2 /v 2 

Here Ii{a \/t 2 — x 2 /v 2 ) is the Bessel function of order 1 : thus —iJiiiz) = 
I\{z). The function is entirely real. The derivation of formula 
(211) is a little troublesome, owing to the discontinuous character of 
the function F: the detailed steps are given in an appendix. 
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The preceding solution depends for its outstanding directness 
and simplicity on the recognition of the infinite integral identity (m), 
into which the integral equation of the problem can be transformed. 
When such identities are known their value in connection with the 
solution of operational equations requires no emphasis. On the 
other hand, we cannot always expect to find such an identity in the 
case of every operational equation; and, particularly in the case of 
such an important case as the transmission equation it would be 
unfortunate to have no alternative mode of solution. Fortunately 
a quite direct series expansion solution is obtainable from the oper- 
ational equation, and this will now be derived. As a matter of con- 
venience we shall restrict the derivation to the voltage formula 

V=e-^ {p+py " T \ ( 203 ) 

As a further matter of mere convenience we shall assume that G = 0, 
so that o = p and (203) becomes 

y =e -Tiip*+2pp (203-a) 

where t = x/v. 

The method holds equally well for the current equation (204) and 
for the general case c^p. 
Write (203-a) as 

y = e -Tp(\+2 P /p) y2 

and expand the exponential factor (l+2p/p) 1/2 by the binomial 
theorem; thus 

a+*/#)«*-i+| +<"{j)'+°"(f)'+ • • • 

so that 

V=e rP.e <".exp( p ~~ p2 p? " ■ ■ • )• 

Now the operational equation 

/ a«Tp 2 a S Tp 3 a t Tp* \ 

can be expanded in inverse powers of p; thus 

P P- P" 

the power series solution of which is 
fA u.^ , ft/ 2 , /3 3 / 3 , 
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It follows at once from the preceding and Theorem VII that 
7(/)=0for Kt 

-.-"(l+ft^+A^ +....) for t>r 

If the coefficients /3i, fii . . . are evaluated, a simple matter of elemen- 
tary algebra, the foregoing expansion in the retarded time t— t will 
be found to agree with the solution (211) when a is put equal to p. 

We shall now discuss the outstanding features of the propagation 
phenomena in the light of equations (210) and (211) for the current 
and voltage. We observe, first, that we have a true finite velocity of 
propagation i»=l/v LC. No matter what the form of impressed 
e.m.f. at the beginning of the line (x = 0), its effect does not reach the 
point * of the line until a time t = x/v has elapsed. Consequently 
v = x/t is the velocity with which the wave is propagated. This is a 
strict consequence of the distributed inductance and capacity of the 
line and depends only on them, since v=\/\/LC. It will be recalled 
that in the case of the cable, where the inductance is ignored, no 
finite velocity of propagation exists. 

The question of velocity of propagation of the wave has been the 
subject of considerable confusion and misinterpretation when dealing 
with the steady-state phenomena. It seems worth while to briefly 
touch on this in passing. 

As has been pointed out in preceding chapters, the symbolic or com- 
plex steady-state formula is gotten from the operational equation by 
replacing the symbol p by ico where i = \/ — 1 and a>/27r is the frequency. 
If this is done in the operational equation (203) for the voltage, the 
symbolic formula is 



V = e — V(»w+p)' o*giat 

If the expression \/(iw+p) 2 — a 2 is separated into its real and im- 
aginary parts we get an expression of the form 

7 =e -a*gto(< p^) f 



WhCre a _ . M+cr 2 -p 2 +V(o; 2 + ( r 2 -p 2 ) 2 +4 t oV 

P V 2co 2 

and 

a = p/f3v. 

X 

Now if we keep the expression t— /3 — constant, that is, if we move 

v 

along the line with velocity dx/dt = v/fi, the phase of the wave will 
remain constant. This is interpreted often as meaning that the 
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velocity of propagation of the wave is v/(3. Now since /3 is greater 
than unity and only approaches unity as the frequency becomes 
indefinitely great, the inference is frequently made that the velocity 
of propagation depends upon and increases to a limiting value v, 
with the frequency. This velocity, however, is not the true velocity 
of propagation, which is always v, but is the velocity of phase propagation 
in the steady -state. This distinction is quite important and failure 
to bear it in mind has led to serious mistakes. 



Returning to equation (211) and (210) we see that after a time 
interval t = x/v has elapsed since the unit e.m.f. was impressed on the 
cable, the voltage at point x suddenly jumps from zero to the 
value e~ px,v while the current correspondingly jumps to the value 

~C - L 



s 



e p-v/r Yhe exponential factor px/v is 



*(£+£)vxMi>lM>l§)-«* 



which will be recognized as the steady-state attenuation factor for high 
frequencies. Similarly -\/C/L is the steady-state admittance of the 
line for high frequencies. The sudden jumps in the current and 
voltage at time t = x/v are called the heads of the current and voltage 
waves. If, instead of a unit e.m.f., a voltage /(/) is impressed on the 
line at time / = 0, the corresponding heads of the waves are f(o)e~ ax 
and \/C/L j{o)e~" x for voltage and current respectively. These 
expressions follow at once from the integral formula 



■^=djy- T)h ^ dT 

=f(o)h(t)+£f'(t-T)h(r)dr. 



The tails of the waves, that is, the parts of the waves subsequent 
to the time t = x/v, are more complicated and will depend on the 
distance x along the line and on the line parameters p and a. The two 
simplest cases are the non-dissipative line, and the distortionless line. 

The ideal non-dissipative line, quite unrealizable in practice, is one 
in which both R and G are zero. In this case p = <r = 0, and formulas 
(210) and (211) become 

1 = for t<x/v, 



[C 

V=0 for t<x/v, 
= 1 for t>x/v. 
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Both current and voltage jump, at time t = x/v, to their steady values. 
If an e.m.f./(/) is impressed on the line at time / = 0, the corresponding 
current and voltage waves are 

1 = for t<x/v, 
= -yj^f(t-x/v)iort>x/v, 

7=0 for t<x/v, 

=f{t-x/v) for t>x/v. 

Consequently the ideal non-dissipative line transmits the waves 
with finite velocity v, without attenuation or distortion. Such a line 
is, of course, the ideal transmission system. 

The non-dissipative line is, of course, purely theoretical and un- 
realizable in practice; the distortionless line is, however, approxi- 
mately realizable, and as the name implies, transmits without distortion 
of wave form. The distortionless line is one in which the line con- 
stants are so related that 

R G a 

If this condition is satisfied, formulas (210) and (211) become 
1 = for t<x/v, 

= J-e~ ax forf>x/v, 

7=0 for t<x/v, 
= e~ ax for t>x/v. 

Furthermore, if the impressed e.m.f. is j(t), the corresponding current 
and voltage waves are: — 

1 = for Kx/v, 
= AQ.e- ax f(t-x/v) for t>x/v, 

V=0fort<x/v, 

= e-" x f(t-x/v) for t>x/v. 

The distortionless line, therefore, transmits the waves without dis- 
tortion of wave form, but attenuates the waves by the factor e~ GX . 
Such a line is an ideal transmission system as regards preservation 
of wave form, but introduces serious attenuation losses. For example, 
if a line has normally negligible leakage, and leakage is introduced 
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to secure the condition R/L = G/C, the line is thereby rendered dis- 
tortionless but the attenuation is doubled. 

One of Heaviside's most important contributions to wire trans- 
mission theory was to point out the properties of the distortionless 
line, its approximately realizable character, and to base on it a correct 
theory of telephonic transmission. 

The character of the wave propagation when the parameters 
p and o- are not restricted to special values, can only be roughly in- 
ferred from inspection of the formulas, and then only when the prop- 
erties of the Bessel function I and Ii have been studied. Fortunately 
these functions have been computed and tabulated for small values 
of the argument, and have simple asymptotic expansions for large 
values. It is therefore a simple matter to compute and graph a 
representative set of curves which show the current and voltage 
waves for various values of p, a and x. For this purpose it is con- 
venient to introduce a change of variables and write : 

r=vt 

a = p/v 

b = <r/v 

whence the formulas for current and voltage become : 

V=e-~+lKf Te - atI f V ^d, (211a) 

Figs. (9) to (18) give a representative set of curves illustrating 
the form of the propagated current and voltage waves for different 
lengths of line, and different values of the line parameters a and b, 
or p and a. 

The curves of Figs. (9) and (10) show the current entering the 

line in response to a unit e.m.f. applied at time / = 0. The line is 

assumed to be non-leaky (6 = 0) and is computed for two different 

values of the parameter a. We see that the current instantly jumps to 

the value y/C/L and then begins to die away, the rate at which Jt dies 

R \C 
away depending on and increasing with the parameter a = — -*l — . 

If we now consider a point x out on the line, the current is zero 
until t = x, at which time it jumps to the value -y/C/L e~ ax . It then 



(210a) 



74 



BELL SYSTEM TECHNICAL JOURNAL 



begins to die away provided x and a are such that ax <2. If, however, 
we are considering a point at which ax>2, the current begins to rise 
instead of fall after the initial jump, and may attain a maximum 
value very large compared with the head before it starts to die away. 
This is shown in the curves of Figs. (11), (12) and (13), also computed 
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Fig. 9 — Current entering line; \— =a=0.O132; G = 0. 
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Fig. 10 — Current entering line; -y- v— =a =0.2645; (7 = 0. 
Multiply orclinates by yjC/L 



for the non-leaky line (6 = 0). From these curves we see that, as 
the length of the line and the parameter a increase, the relative mag- 
nitude of the tail, as compared with the head of the wave, increases. 
Finally when the line becomes very long, the head of the wave be- 
comes negligibly small, and the wave, except in the neighborhood of 
its head, becomes very close to that of the corresponding non-inductive 
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cable. This is shown in curves (13) and (14), for the line and the 
corresponding cable, which are plotted to the same time scale and 
ordinate scale to facilitate comparison. Curve (15) shows the effect 
of leakage in eliminating the tail. This line is not quite distortionless 
but nearly so. 
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Fig. 11— Propagated current inline; .v = 200; ■s-'Vt" a = °- 01 32; G = 0. 
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An interesting feature of both current and voltage waves is that 
when a sufficient time has elapsed after the arrival of the head of the 
wave, the waves become closer and closer to the wave of the cor- 
responding non-inductive cable; that is, to the cable having the same 
R,C and G. Consequently the inductance plays no part in the 
subsidence of the waves to their final values. 
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Fig. 13 — Propagated current in li^v, _ , T 
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Fig. 14 — Propagated current in cable. Multiply ordinates by 2/22*. 
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Curves (16), (17) and (18) illustrate the voltage wave for several 
conditions. After the arrival of the head, the wave slowly builds 
up to its final value. Curve (18) represents the case where the line 
is very nearly distortionless, showing how completely the distorting 
tail of the wave is eliminated. 
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Fig. 16 — Propagated voltage in line; — "\l y * = a* = 3 ; ^ = 0. 

So far we have confined attention to the current and voltage waves 
in response to a unit e.m.f. applied at time / = to the line terminals. 
Of much greater technical importance is the question of the waves 
in response to a sinusoidal e.m.f. suddenly applied to the line termi- 
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nals. In order to investigate this important problem it is convenient 
to divide the expressions for the current and voltage waves as given 
by equations (210-a) and (211-a) into two components. We write 
for t>x, 

"c 



1 = 



L 



e~ a * +/(/), 



V=e-"*+W(l), 



(210-b) 
(211-b) 



where, by definition, J(t) and W{t)' are the differences between the 
total waves and their heads. The advantage of analyzing the waves 
into these components is that the distortion of the waves is due to 
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/(/) and W(t) respectively, while the first component of (210-b) and 
(211-b) introduce merely a delay. Thus, if the e.m.f. impressed 
at time / = is/(/), the corresponding waves for />.v/y or t~2:x, are 



= ^~e-°\f(t-x/v) 

Jx/v 



V = e-" x f(t-x/v) 



+ I f{t-h)W'{h)dh, 

Jx/v 



(212) 



(213) 



where J'(t) =^JQ) and W'{t) =~ W{t). 

The integrals of (212) and (213) can be computed and analyzed 
in precisely the same way as discussed in connection with the non- 
inductive cable problem, and are of very much the same character 
as the alternating current waves of the cable. In the total waves, 
however, as given by (212) arid (213), a very essential difference is 
introduced by the absence of the first terms, which represent undis- 
torted waves propagated with velocity v. Thus, if the impressed 
e.m.f. is sin ut, (212) and (213) become 



/= -%! j-e " x sin <x)(t — x/v) 

4. / sin u(t—ti)J'(ti)dti, fort>x/v 

Jx/v 



(214) 



V=e~ ax sin u(t—x/v) 



(215) 
sin w{t—h)W'(h)dh, \ovt>x/v. 

x/v 



4- / sin u(t-t 1 )W'(t 1 )dt 1 , for t>x/v. 

tJx/v 



Now the first terms of (214) and (215) are simply the usual steady- 
state expressions for the current and voltage waves when the fre- 
quency is sufficiently high to make the steady-state attenuation 
constant equal to a and the phase velocity equal to v. Furthermore 
the integral terms become smaller and smaller as the applied fre- 
quency w/27r is increased. It follows, therefore, that for high fre- 
quencies the waves assume substantially their final steady value 
at time t = x/v, and that the tails of the waves, or the transient 
distortion, becomes negligible. This is a consequence entirely of the 
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presence of inductance in the line, and shows its extreme importance 
in the propagation of alternating waves and the reduction of transient 
distortion. 

It should be pointed out, however, that if the line is very long and 
the attenuation is very high, the integral terms of (214) and (215) 
are not negligible unless the applied frequency is correspondingly very 
high. For example, on a long submarine cable, the a.c. attenuation 
is so large that the first terms of (214) and (215) are very small, and 
J(t) is very large compared with \/C/L e~ ax . Consequently here 
there is very serious transient distortion and alternating currents 
are therefore not adapted for submarine telegraph signalling. 

This discussion may possibly be made a little clearer, without 
detailed analysis, if we recall the discussion of alternating current 
propagation in the non-inductive cable of the preceding chapter. 
From that analysis it follows that, when the applied frequency tc/2ir 
is sufficiently high, the integral term of (214) becomes approximately 

and the complete current wave is 

■\jj-e~ ax sin co(t-x/v) +—/'(*) (216) 

and similarly the voltage wave is 

e~ ax sin to(t-x/v) + - W'(t). (217) 

CO 

Now if the total attenuation ax is large the last terms of (216) and 
(217), before they ultimately die away, may become very large com- 
pared with the first terms, which represent the ultimate steady-state. 

Appendix to Chapter VII. Derivation of Formula (211) 

The only troublesome question involved in deriving (211) from 
(207) and (209) is that we have to differentiate with respect to x, in 
accordance with (207), the discontinuous function F(t). To accom- 
plish this we write (209) in the form 

F(t) =<j>(t-x/v)e-P'I (trVt 2 -x 2 /v 2 ) (209-a) 

where <j>(t) is defined as a function which is zero for t<x/v and unity 
for t>x/v. Clearly this is equivalent to (209) and permits us to deal 
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with F(t) as a continuous function. Now, in accordance with (207), 
perform the operation of differentiation upon (209-a) : we get 

- V ^L = %<j,(t-x/v)e-P<U<r\/¥=¥M 
ox ot 

-v4>(t-x/v) J- e-vioiffVP-xW. 

The first expression follows from the fact that 

We observe also that ^ <f>(t-x/v) =0 except at t = x/v, when it is in- 

ot 

finite. We also observe that, for t>x/v, 

£^(t-x/v)dt = l 

and that the whole contribution to the integral occurs at t = x/v. 
With these points clearly in mind, the expression 

Jo dx 
reduces to (211) without difficulty. 

CHAPTER VIII 

Propagation of Current and Voltage in Artificial 
Lines and Wave Filters 

The artificial line here considered is a periodic structure, com- 
posed of a series of sections connected in tandem, each section con- 
sisting of a lumped impedance Z\ in series with the line, and a lumped 



Fig. 19 



impedance s 2 in shunt across the line. In the artificial line which 
we shall consider it will be assumed that the voltage is applied at 
the middle of the initial or zeroth section, as shown in Fig. 19. This 
termination is chosen because of its practical importance, and be- 
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cause also of the fact that the mathematical analysis is simplified 
thereby. Furthermore any other termination can be regarded and 
dealt with as an additional terminal impedance, so there is no essen- 
tial loss of generality involved. 

A study of the properties of the artificial line is of practical im- 
portance for several reasons : 

1. The artificial line is often used as a model of an actual trans- 
mission line and it is therefore of importance to determine theoretically 
the degree of correspondence between the two. 

2. The solution for the corresponding transmission line with con- 
tinuously distributed constants is derivable from the solution for the 
artificial line by keeping the total inductance, resistance, capacity 
and leakage constant or finite, and letting the number of sections 
approach infinity. 

3. The artificial line is very closely related, in its properties and 
performance, to the periodically loaded line, and its solution is, to 
a first approximation, a working solution for the loaded line. 

4. The structure is of great importance in its own right, and when 
the impedance elements are properly chosen, constitutes a "wave 
filter." 

We shall now derive the operational and symbolic equations which 
formulate the propagation phenomena in the artificial line. Let I„ 
denote the mesh current in the nth section of the line; I,,-i the mesh 
current in the (it — l) th section, etc. Now write down the expression 
for the voltage drop in the n th section; in accordance with Kirch- 
hoff's law we get : 

( 2l +2s 2 )/ / ,-s 2 (/„- 1 + /„ +1 )=0 (218) 

where, of course, the impedances have the usual significance. 

Now this is a difference equation, as distinguished from a differ- 
ential equation, but the method of solution is essentially the same. 
We assume a solution of the form 

I„=Ae-»v+Be'* (219) 

where A, B and T are independent of n, and substitute in (218). 
After some simple rearrangements we get 

\ ( 2l +2s 2 )-2z 2 cosh Y\ ■ [Ae-T+Be*?] =0. (220) 

Equation (218) is clearly satisfied by the assumed form of solution, 
and furthermore leaves the constants A and B arbitrary and at our 
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disposal to satisfy any boundary conditions, provided V is so chosen 
that 

cosli I = — -= , „ v 

2s 2 (221) 

where p = Zi/42 2 . 

Now by reference to equation (219) it is easily seen that T is the 
propagation constant of the artificial line, precisely analogous to the 
propagation constant y of the smooth line. In terms of the im- 
pedances Zi and 2 2 , the propagation constant of the artificial line is 
determined by (221). This equation may either be regarded as an 
operational equation or a symbolic equation, depending on whether 
the impedances are expressed in terms of the operator p or in terms 
of iu, where co is 27r times the frequency. 

Now suppose in (221) we write e v = x; the equation becomes 

.v+l/.v = 2(l + 2 P ) 
and solving for x we get 

* = < J r = (l + 2p) + V(l+2 P , 2 -l 

=(ViTp+x / p) 2 =(Vi+p-Vp) 2 

which is an explicit formula for T. 

Now return to equation (219) and let us assume that the line is 
either infinitely long, or, what amounts to the same thing, that it is 
closed by an impedance which suppresses the reflected wave. We 
assume also that a voltage V is impressed at mid-series position of the 
zero"' section (»=;). Equation (219) becomes 

I a =Ae~ n ? 

and the currents in the zero"' and 1*' sections are 

l =A, h=Ae-r. 

Now, by direct application of Kirchhoff's law to the zero'' section, 
we have 

whence 

A)hz i +z,(l-e-n\ = V . (223) 

But 

_ r =A=~v Ol 

I = W-r 
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where K is the characteristic impedance of the artificial line (at mid- 
series position). Hence by (223) and (222) 

11 1 

K z 2 (l-e _r )+2p 

(224) 
= _1 3_ = l_ 1_ 

2z 2 \/p-\- P - \/ziZi s/l +p 
By aid of the preceding the direct current wave can be written as 

j = V J_ tVl+P-Vpf ^ (225) 

Vzi22 Vl+P 

This formula is not so physically suggestive as its equivalent 

but is useful when we come to the solution of the operational equation. 
Before proceeding with the operational equation, and the investi- 
gation of transient phenomena in artificial lines, it will be of interest 
to deduce from the foregoing the unique and remarkable properties 
of wave filters in the steady state. For this purpose we return to 
equation (221) 

cosh r = l+2p. 

Now suppose that the series impedance Z\ is an inductance L and 
the shunt impedance z 2 a capacity C, so that, symbolically, 

. T 1 co 2 LC 

tut 4 

and 

cosh r = l-ico 2 LC. (226) 

Now let us write T = i9, where i = \/—l\ the preceding equation 
becomes 

cos = 1-| «?LC (227) 

and the ratio of currents in adjacent sections is e~ t9 . Consequently if 
6 is a real quantity the ratio of the absolute values of the currents in ad- 
jacent sections is unity, and the current is propagated without attenuation. 
Inspection of equation (227) shows that 6 is real provided the right 
hand side lines between +1 and —1: or that co lies between and 
2/\/'LC. Consequently this type of artificial line transmits, in the 
steady state, sinusoidal currents of all frequencies from zero to 
l/irs/LC without attenuation. It is known as the low-pass filter. 
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If we invert the structure, that is, make the series impedance Z\ 
a capacity C and the shunt impedance z 2 an inductance L, so that 

1 . T 1 

we get, corresponding to (226) and (227), 

COshr = 1 -2^C' (228) 

cosd = 1 -2^LC- < 228a > 

This type of artificial line transmits without attenuation currents 
of all frequencies for which the right hand side of (228-a) lies between 
+ 1 and — 1; that is, all frequencies from infinity to a lower limiting 
frequency \/4.v\/LC, while it attenuates all frequencies below this 
range. It is known, on this account, as the high-pass filter. 

It is possible by using more complicated impedances to design 
filters which transmit a series of bands of frequencies. We cannot, 
however, go into the complicated theory of wave filters here, which 
has been covered in a series of important papers. One point should 
be noted, however: transmission without attenuation implies that 
the impedance elements are non-dissipative. Actually, of course, 
all the elements introduce some loss, so that in practice the filter 
attenuates all frequencies. Careful design, however, keeps the 
attenuation very low in the transmission bands. 

We shall now derive the indicial admittance formulas for some 
representative types of artificial lines and wave filters from the oper- 
ational formula 

An = .*_ h/l+P + V7r 2 "- (229) 

V(1+p)ziZ2 

This equation follows directly from (225) on putting V = l. 

We start with the so-called low-pass filter on account of its sim- 
plicity and also its great importance in technical applications. This 
type of filter consists of series inductance L and shunt capacity C. 
The general case which includes series resistance R and shunt leakage 
G has been worked out (see Transient Oscillations, Trans. A. I. E. E., 
1919). The solution is, however, extremely complicated and will not 
be dealt with here. We shall, instead, consider the important and 
illuminating case where the series and shunt losses are so related 
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as to make the circuit quasi-distortionless. We therefore take, 
operationally, 

z l =pL+R = L{p+\) 

(230) 
l/zi=pC+G=C(p+\) 

where X=R/L = G/C. 
We then have 

2l2 2 = L/C, 

*/— *W+M (23]) 

P-^+X)». 

Now by reference to formula (229) we see that >1„ is a function of 
(p+\)', thus 



^ H = ^i = f 1+ M 



Z w (£+X) V p/(p + X)Z„( P + X)- 
Now write 

4-= -L- 
" z„(/>)' 

It follows at once from reference to theorem VII that 

A n =(l+\rdt}A°e~ U (232) 

so that the problem is reduced to the solution of the operational 
equation for A",. Writing o) c = 2/\ / LC, we have 

£ Vi+\p-mI Vi+5753*+>/*] "" 

(233) 

|c ^ rVj> 2 +q>c 2 -^ -| 2 " 

Now refer to formula (w) of the table of integrals; writing y/L/C=k, 
we see by Theorem Y that 



. „ 1 ('uj r / s , 



(234) 



where Jm(j) is the Bessel function of order 2n and argument r. We 
note also that this is the indicial admittance of the non-dissipative 
low-pass wave filter; that is, the current in the w th section in response 
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to a unit e.m.f. applied to the initial section (w = 0). From (232) 
and (234) it follows at once that 



. ft fO) c T 

r / dre' XT / JUr^dr,. 
ija Jo 



Integrating the second member by parts and noting that A°,{o)=Q, 
this reduces to 

1 pu e t \ 



A 



1 C "' 1 — 
n= k \ e ^/ J,„(T)dT (235) 



which is the indicial admittance formula for the quasi-distortionless 
low-pass filter, or artificial line. 

Before discussing these formulas, it is of interest to derive the 
formula for A", by power series expansion. Formula (233) can be 
written 



" k \pJ vi + W^li+Vi+k/M 2 "' 

This can be expanded in a series in inverse powers of p; thus 
1 ( /u)c\ 2 " +{ 2ii + 2 /co cN 2 '+ 3 



A " = k2 2 "\ \p) 2 2 1! \p) 



(2rc + 3) (2rc+4) (ok\ 2 " +b _ 
+ 2'2! \p) 

Replacing l/p" by t"/n\ in accordance with the Heaviside Rule we get 

A „_l_ f _1 (oht\ 2 "+ ' _ _ 2w+2 fuj \ 2 "+ 3 

"~ k }(2«+l)!\ 2 ) l!(2« + 3)!\2 ) 

(2n+3 )(2n+4) My»+* 1 (235 „ a) 

^ 2!(2ra + 5)! \2) j 

This can be recognized as the power series expansion of (234). 

The artificial cable is also of interest and practical importance. 
In this structure the series impedance is a resistance R and the shunt 
impedance is a capacity C, so that 

Z!=R, l/zo=pC, 
z iZ2 =R/pC, zJz^pRC, (236) 

P = pRC/4. 
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Now let us return to formula (229), and expand in inverse powers of 
p : we get 

A __ 1 ( 1 2n+2 1 (2w+3)(2w+4) 1 _ I 

n= 2W^iJ~? 22 1! P" +l 2 <2! P" +2 ••••f ( ^ 7) 

Now since v PZ1Z2 = — , we have 



A __?_JY 2 \» 2w+2 / 2 \»+i 
n ~ 2"R\ \RCp) 2.1!" V-RcW 



(2«+3)(2»+4) / 2 \"+ 2 
+ 2 2 2! Vi?C^J 

Replacing 1/p" by r/w! we get finally 

A 2 fl/2/\» (2«+2) (2t\"+ l 

An 2»R\n\\RC) 2.V.{n + l)\\RC) 

( 2w+3)(2w+4) /2/ \ »+ 2 _ 
+ 2 2 .2!(w+2)! \RC) ' 



(238) 



For large values of n and t this series is difficult to compute or in- 
terpret. It can, however, be recognized as the series expansion of the 
function 

'-!'"**(£) (239) 

where I„(2t/RC) is the Bessel function /„ of order n and argument 
(2t/RC). This solution, it may be remarked, can be derived directly 
by a modification of the integral formula («). 

It is beyond the scope of this paper to consider other types of 
artificial lines and wave filters; for a fairly extensive discussion the 
reader is referred to "Transient Oscillations in Electric Wave-Filters," 
B. S. T. J., July, 1923. The low-pass wave filter, however, both in 
its own right and on account of its close relation to the periodically 
loaded line, deserves further discussion. 

For the non-dissipative low-pass wave filter, we have 

A° n = \j" ct J 2n {T)dr (234) 

while for the quasi-distortionless low-pass wave filter 

A n =\j^ ct e-^J in {r)dr (235) 

where p. = \/co c = R/Lco c = R/2vL. 
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Computation and analysis of these formulas involve an elementary 
knowledge of Bessel functions. The properties necessary for our 
purposes are briefly discussed in an appendix to this chapter. 

The indicial admittances for the non-dissipative low-pass filter, 
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Fig. 20 — Low pass wave filter. Indicial admittance of initial section (»=0). 
Multiply ordinates by VC/L 
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Fig. 21 — Low pass wave filter. Indicial admittance of third section (» = 2). 
Multiply ordinates by VC/L 
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Fig. 22 — Low pass wave filter. Indicial admitt ance of fifth section (w=4). 
Multiply ordinates by VC/£ 

that is, the current in response to a steady unit e.m.f. applied at 
time / = 0, are shown in the curves of Figs. 20, 21 and 22, for the 
initial or zeroth, the 3rd and the 5th sections, respectively. These 
curves together with the exact and approximate formulas given 



on 



BELL SYSTEM TECHNICAL JOURNAL 



above are sufficient to give a reasonably comprehensive idea of the 
general character of these oscillations and their dependence on the 
number of sections and (he constants of the filter. 

It will be observed that the current is small until a time approxi- 
mately equal to 2n/u>, - = ny/h\C% has elapsed after the voltage is 
applied. Consequently the low-pass filter behaves as though cur- 
rents were transmitted with a finite velocity of propagation a>,/2 = 
1/vZ/iCo sections per second. This velocity is, however, only ap- 
parent or virtual since in every section the currents are actually 
finite for all values of time>0. 

After time t = n\/ L\C* has elapsed the current oscillates about the 
value 1/k with increasing frequency and diminishing amplitude. 
The amplitude of these oscillations is approximately 



1/k 



Vl-(2n/co c t) 



*l 



2 

7TOV/ 



and their instantaneous frequency (measured by intervals between 

zeros) 

(tic , 

2^Vl-(2»/«,/) 2 . 



The oscillations are therefore ultimately of cut-off or critical frequency 
(a c /2ir in all sections, but this frequency is approached more and more 
slowly as the number of filter sections is increased. 

Figs. 23, 24, 25, give the indicial admittance in the 100th, 500th 
and 1000th section of the low-pass wave filter. The filter itself seldom 
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Fig. 23 — Low pass wave filter. Indicial admittance of 100th section (// =99). 
Multiply ordinates by yJi/L 
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embodies more than 5 sections. The case of a large number of sec- 
tions is of interest, however, because it represents a first approxi- 
mation to the periodically loaded line. While the non-dissipative 




1000 1020 1040 1060 1030 1100 1120 1140 1100 

Fig. 24 — Low pass wave filter. Inclicial admittance of 500th section (n=499). 
Multiply ordinates by VC/L 




Fig. 25 — Low pass wave filter. Indicia! admittance of 1000th section (« = 999). 
Multiply ordinates by VC/Z, 



line is ideal and unrealizable, its study is of practical importance 
because in this type of line the effect of the discontinuous character 
of the loading of the periodically loaded line is isolated and exhibited 
in the clearest possible manner. 
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The dotted curves represent the current in the corresponding 
smooth line. For the smooth line, the current, as we have seen, is 
discontinuous, being identically zero for a_time vt = n and having an 
instantaneous jump to its final value y/C/L at vt = n. The current 
in the artificial or periodically loaded line differs from that in the cor- 
responding smooth line in three important respects: (1) the absence 
of the abrupt discontinuous wave front, (2) the presence of super- 
posed oscillations, and (3) the absence of a true finite velocity of propaga- 
tion. It will be observed, however, that the current in any section 
is negligibly small or even sensibly zero until vt = n, so that the current 
is propagated with a virtual velocity l/\/LC per section. The pres- 
ence of a well marked wave front is also evident although this is not 
abrupt, as in the smooth line. The effective slope of the wave front 
becomes smaller as the current wave travels out on the line, decreasing 
noticeably as the number of sections is increased. When the number of 
sections becomes large, however, the decrease in the slope is not rapid, 
being in the 500 th section about 60 per cent, of that in the 100 th section. 

The superposed oscillations are of interest. These are initially 
of a frequency depending upon and decreasing with the number of 
sections, n, but in all sections ultimately attaining the frequency 

1_ = v_ 

■K\/LC 7T 

which is the critical or cut-off frequency of the line, above which 
steady-state currents are attenuated during transmission and below 
which they are unattenuated. When vt is large com pared with n 
the amplitude of these oscillations becomes s/l/rvt so that they 
ultimately die away and the current approaches the value V 'C/L 
for all sections. The current in the loaded line is thus asymptotic 
to the current in the corresponding smooth line and oscillates about 
it with diminishing amplitude and increasing frequency. 

Since the abscissas of these curves represent values of 2vt = 2t/\/LC, 
and the ordinates are to be multiplied by s/CfL to translate into 
actual values, the curves are of universal application for all values 
of the constants L and C. 

The investigation of the building-up of alternating currents in wave 
filters and loaded lines is very important. It depends for the non- 
dissipative case on the properties of the definite integrals 

sin WT J„(t)(It, 

XU c t 
COS WT Jn(j)dr, 
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where w = a/<a c and o) = 2tt times the applied frequency. The mathe- 
matical discussion is, however, quite complicated and will not be 
entered into here. The reader, who wishes to follow this further, 
is referred to Transient Oscillations, Trans. A. I. E. E., 1919 and 
Transient Oscillations in Electric Wave Filters, B. S. T. J., July, 1923. 

Appendix to Chapter VIII. Note on Besel Functions 

The Bessel Functions of the first kind, J„(x) and I nix), are defined, 
when n is zero or a positive integer, by the absolutely convergent 
series 



Mx)= 2^n-\\ 



1- 



x 2 



.%•> 



2(2«+2) ' 2.4(2»+2)(2»+4) 
x* 



2.4.6.(2»+2)(2»+4)(2»+6) 



/„(*) = 



x" 
2".n\ 



X 1 



r^ + TTTT, 



.%•> 



+ 



2(2»+2) ' 2.4(2«+2)(2«+4) 

x 6 
2.4.6.(2«+2)(2»+4)(2»+6) "* 



In the following discussion of the properties of these functions it will 
be assumed that the argument x is a pure real quantity. 

For large values of the argument (x large compared with n), the 
behavior of the functions is shown by the asymptotic expansions: — 

4w 2 -l . (4» 2 -l)(4« 2 -9) 



In(x)=-^-=. 

■\Z2ttx 



1- 



l!(8s) ' 2! (8s) 2 
(4w 2 -l)(4w 2 -9)(4« ! -25) 



/ 



n(x) = .J ] P„ COS (X — 



3!(8x) 3 
2w+l 



+ . .. 



if) — Qn sin (x - A 7r) 



where 



4 "' *-— v~ 4 

(4« a -l)(4n 8 -9) (4tt 2 -l)(4tt 2 -9)(4rc 2 -25)(4rc 2 -49) 
" 2!(8*> 2 " 41 (Rx)* •'■' 



2! (8*) 2 ' 4! (Sx)* 

4w 2 -l (4rc 2 -l)(4w 2 -9)(4rc 2 -25) 
8x 3!(8x) 3 + • • • • 

We thus see that I„ increases indefinitely and behaves ultimately as 

g* 

V2lrx' 
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The function J n {x), however, is oscillatory and ultimately behaves as 



I 2 / 2w+l \ 



For all orders of » 

J n (x)dx = l. 



£ 



The properties of / M (x) may be described qualitatively as follows: — 

When the argument is less than the order (0^x<w) the function 
is very small and positive, and is initially zero (except when n = 0). 
In the neighborhood of x = n, the function begins to build up and 
reaches a maximum a little beyond the point x = n. Thereafter the 
function oscillates with increasing frequency and diminishing ampli- 
tude, and ultimately behaves as 



>l 



/ 2»+l \ 

V—T-') 



2 / 2w+l 
— cos I 

WX 



When « = 0, the initial value is unity, but the subsequent behavior of 
the function is as described above. 

A more precise description of the function is gotten from the follow- 
ing approximate formulas. 

Jn (x) = B„ (x) cos fi„ (x) , for x > n 

where 

Bn{x)= X \\ l 



■kx f 1 m 2 , 3 m 

*■ ~ !5T "T" o 



x 2 ' 2 .x- 4 (l-ra 2 /x 2 



/ x f ll w2 ■ W . _./OT\ OT 2 1 "I 



2»+l 

: — 7T, 



tf a (x) = ~Sl n {x), 



=^ 



w 2 , 3ffl s 



and 



x 2 ' 2 x 4 (l-m 2 /x 2 ) 2 ' 

W 2 =K 2_ 1 /4_ 



This approximate formula is valid only where x>n, its accuracy 
increasing with x and with n. For all orders of n it is quite accurate 
beyond the first zero of the function. 
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The "instantaneous frequency" of oscillation is approximately 



1 n't~\ * 1 w2 . 3 m ' 1 

27r"" ( - Vj "27r\ i * 2+ 2 ** (1-m 2 /* 2 ) 2 ' 

By this it is meant that at any point x(x>n) the interval between 
successive zeros is approximately ir/2'(x). Otherwise stated, in 
the neighborhood of any point x, the function behaves like a sinusoid 
of amplitude B„(x) and frequency o/27r where o> = Q„(x). 

The following approximate formulas, while not sufficiently precise 
for the purposes of accurate computation except for quite large values 
of x, clearly exhibit the character of the functions for values of the 
argument x>n, and of the order n>2. 

J„ (x) = h„ \l —'cos (q„x — d„), 

l~2~ 

J' n (-"v) = —qJh, \ — sin (q„x — 6„), 

\ 7T.V 

/ J„(.r)r/.v = l + — a — sin (q„x — d„), 
Jo q virx 

where 

hn -\l-n 2 /xV + 4x 2 ' 



q„ = \ / l-n 2 /x 1 . 
„ 2;/ + l 



t — w sin '(m/.v). 



